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Abstract 



in . 

We use the BLM method to show that perturbatively-calculable observables 
in QCD, including the annihilation ratio R e + e -, the heavy quark potential, 
^p Hi and radiative corrections to structure function sum rules can be related to each 

other without renormalization scale or scheme ambiguity. The commensurate 
scale relations connecting the effective charges for observables A and B have 
the form o^Qa) = ob(Q.b) (l + r A / B ^- + •••), where the coefficient r^j B 
is independent of the number of flavors / contributing to coupling constant 
renormalization. The ratio of scales Qa/Qb is unique at leading order and 
guarantees that the observables A and B pass through new quark thresholds at 
the same physical scale. We also show that the commensurate scales satisfy 
the renormalization group transitivity rule which ensures that predictions 
in PQCD are independent of the choice of an intermediate renormalization 
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scheme C. In particular, scale-fixed predictions can be made without reference 
to theoretically-constructed renormalization schemes such as MS. QCD can 
thus be tested in a new and precise way by checking that the observables 
track both in their relative normalization and in their commensurate scale 
dependence. 

The generalization of the BLM procedure to higher order assigns a differ- 
ent renormalization scale for each order in the perturbative series. The scales 
are determined by a systematic resummation of running coupling constant 
effects. The application of this procedure to relate known physical observ- 
ables in QCD gives surprisingly simple results. In particular, we find that 
up to light-by-light type corrections, all terms involving £3, £5 and tt 2 in the 
relation between the annihilation ratio R e + e - and the Bjorken sum rule for 
polarized electroproduction are automatically absorbed into the renormaliza- 
tion scales. The final series has simple coefficients which are independent of 
color: a gi (Q) = a R (Q*) - a 2 R (Q**) + a R {Q***), where a = (3C f /4tt)q. The 
coefficients in the commensurate scale relation can be identified with those 

obtained in conformally-invariant gauge theory. 
ll.10.Gh, 11.15.Bt, 12.38.Bx, 13.65.+i 



2 



Typeset using REVTgX 



I. INTRODUCTION 



One of the most serious difficulties preventing precise tests of QCD is the scale ambiguity 
of its perturbative predictions. Consider a measurable quantity such as p = R e + e - (s) — 3Se^. 
The PQCD prediction is of the form 



p = r a s (p) 



i + n(ji) + r 2 (p) -^p + 

7T 7T Z 



Here a s (p) = g 2 /4:7r is the renormalized coupling defined in a specific renormalization scheme 
such as MS, and p is a particular choice of renormalization scale. Since p is a physical quan- 
tity, its value must be independent of the choice of p as well as the choice of renormalization 
scheme. Nevertheless, since we only have truncated PQCD predictions to a given order af , 
the predictions do depend on p. In the specific case of R e + e -, where we have predictions 
through order a^, the sensitivity to p has been shown to be less than 10% over a 
large range of Inp 0. However, in the case of the hadronic beauty production cross sec- 
tion (da/dp^,)(pp — > B + X), which has been computed to next-to- leading order in a s , the 
prediction [|J for the normalization of the heavy quark distribution at hadron colliders 



ranges over a factor of 4 if one chooses one "physical value" such as p = \ J m 2 B + p\ rather 



than an equally well motivated choice p = ym 2 3 +p\. 

There is, in fact, no consensus on how to estimate the theoretical error due to the scale 
ambiguity, what constitutes a reasonable range of physical values, or indeed how to identify 
what the central value should be. Even worse, if we consider the renormalization scale p 
as totally arbitrary, the next-to-leading coefficient r\{p) in the perturbative expansion can 
take on the value zero or any other value. Thus it is difficult to assess the convergence of the 
truncated series, and finite-order analyses cannot be meaningfully compared to experiment. 

The p dependence of the truncated prediction p N is often used as a guide to assess the 
accuracy of the perturbative prediction, since this dependence reflects the presence of the 
uncalculated terms. However, the scale dependence of pn only reflects one aspect of the total 
series. This point has also been recently emphasized by Maxwell et al. H. For example, 
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consider the orthopositronium J PC = 1 decay rate computed in quantum electrodynamics: 
T = T [1 — 10.3 (a/n) + •••]. The large next-to-leading coefficient, r 1 = 10.3 shows that 
there is important new physics beyond Born approximation. The magnitude of the higher 
order terms in the decay rate is not related to the renormalization scale since the QED 
coupling a does not run appreciably at the momentum transfers associated with positronium 
decay. 

Thus we have a difficult dilemma: If we take fi as an unset parameter in PQCD pre- 
dictions, then we have no reliable way to assess the accuracy of the truncated series or the 
parameters extracted from comparison with experiment. If we guess a value for /i and its 
range, then we are left with a prediction without an objective guide to its theoretical pre- 
cision. The problem of the scale ambiguity is compounded in multi-scale problems where 
several plausible physical scales enter. 

In fact three quite distinct methods to set the renormalization scale in PQCD have been 
proposed in the literature: 

1. Fastest Apparent Convergence (FAC) [H. This method chooses the renormalization 
scale p so that the next-to-leading order coefficient vanishes: ri(p) = 0. 

2. The Principle of Minimum Sensitivity (PMS) 0. In this procedure, one argues that the 
best scale is the one that minimizes the scale dependence of the truncated prediction 
Pn, since that is a characteristic property of the entire series. Thus in this method one 
chooses /i at the stationary point dpN/dp = 0. 

3. Brodsky-Lepage-Mackenzie (BLM) 0. In the BLM scale-fixing method, the scale 
is chosen such that the coefficients are independent of the number of quark fla- 
vors / renormalizing the gluon propagators. In leading order, one chooses the 
scale so that / does not appear in the next-to- leading order coefficient. That is, if 
ri(/x) = riotju) +rn(yu)/, where rio(/-t) and ru(p) are / independent, then one chooses 
the scale p given by the condition rn(p) = 0. This prescription ensures that, as in 
quantum electrodynamics, vacuum polarization contributions due to fermion pairs are 
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all incorporated into the coupling rather than the coefficients. In the case of 
non-Abelian theory, the BLM method automatically resums the corresponding gluon 
as well as quark vacuum polarization contributions since the coupling a s is a function 
of A, a 11- §/. 



FIG. 1. The scale [ij^fs according to the BLM (dashed-dotted), PMS (dashed), FAC (full) and 
^Jy (dotted) procedures for the three-jet rate in e + e~ annihilation, as computed by Kramer and 
Lampe || . Notice the strikingly different behavior of the BLM scale from the PMS and FAC scales 
at low y. In particular, the latter two methods predict increasing values of fi as the jet invariant 
mass M. < \/\ys) decreases. 

These scale-setting methods can give strikingly different results in practical applications, 
For example, Kramer and Lampe have analyzed || the application of the FAC, PMS, and 
BLM methods for the prediction of jet production fractions in e + e~ annihilation in PQCD. 
Jets are defined by clustering particles with invariant mass less than y/ys, where y is the 
resolution parameter and \fs is the total center-of-mass energy. Physically, one expects the 
renormalization scale fi to reflect the invariant mass of the jets, that is, \i should be of order 
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y/ys. For example, in the analogous problem in QED, the maximum virtuality of the photon 
jet which sets the argument of the running coupling a(Q) cannot be larger than ^fys. Thus 
one expects fi to decrease as the resolution parameter y — > 0. However, the scales chosen 
by the FAC and PMS methods do not reproduce this behavior (see Fig. 1): The predicted 
scales fipMs(y) and fiFAc(y) rise without bound at small values for the jet fraction y. On the 
other hand, the BLM scale has the correct physical behavior as y — > 0. Since the argument 
of the running coupling becomes small using the BLM method, standard QCD perturbation 
theory in a s [fiBLM{y)] will not be convergent in the low y domain ||. In contrast, the scales 
chosen by PMS and FAC give no sign that the perturbative results break down in the soft 
region. 

In this paper we shall use the BLM method to show that all perturbatively-calculable 
observables in QCD, including the annihilation ratio R e + e -(Q 2 ), the heavy quark potential, 
and the radiative corrections to the Bjorken sum rule can be related to each other at fixed 
relative scales. The "commensurate scale relation" for observables A and B in terms of their 
effective charges has the form 

<xa{Qa) =u b (Qb) (l+r j4/B ^ + ---) • (1.2) 

The ratio of the scales Xa/b = Qa/Qb is chosen so that the coefficient t^/b is independent of 
the number of flavors / contributing to coupling constant renormalization. This guarantees 
that the effective charges for the observables A and B pass through new quark thresholds 
at the same physical scale. We shall show that the value of \a/b is unique at leading order, 



and that the relative scales satisfy the transitivity rule [10 



^a/b — ^a/c ^c/b ■ (1-3) 



This is equivalent to the group property defined by Peterman and Stiickelberg [|XT|| which 
ensures that predictions in PQCD are independent of the choice of an intermediate renormal- 
ization scheme C . In particular, scale- fixed predictions can be made without reference to 
theoretically-constructed renormalization schemes such as MS; QCD can thus be tested by 
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checking that the observables track both in their relative normalization and commensurate 
scale dependence. O 



II. COMMENSURATE SCALE RELATIONS 



It is interesting that the task of setting the renormalization scale has never been consid- 
ered a problem or ambiguity in perturbative QED. For example, the leading-order parallel- 
helicity amplitude electron-electron scattering has the form 



Here a{Q) = a(Q )/(l — U[Q 2 ,Q 2 l ,a(Qo)}) is the QED running coupling which sums all 
vacuum polarization insertions II into the renormalized photon propagator. The value a(0) 
is conventionally normalized by Coulomb scattering at t = —Q 2 = 0. Notice that both 
physical scales t and u appear as arguments of the running coupling constant in the various 
terms contributing to the scattering cross section; if one chooses any other scale for the 
running coupling, in either the direct or crossed graph amplitude, then one generates a 
spurious geometric series in / (a / V) £n(—t / /i 2 ) or / (a/n) £n(—u/ fi 2 ), where / represents 
the number of fermions contributing to the vacuum polarization of the photon propagator. 

In general, the "skeleton" expansion of Feynman amplitudes in QED guarantees that 
all dependence of an observable on the variable / is summed into the running coupling 
constant; the coefficients in QED perturbation series are thus always /-independent once 
the proper scale in a has been set. Note that the variable / is defined to count only vacuum 
polarization insertions, not light-by-light loops, since such contributions do not contribute 
to the coupling constant renormalization in QED. 

The use of the running coupling constant a(Q) in QED allows one to sum in closed form 
all proper and improper vacuum polarization insertions to all orders, thus going well beyond 
ordinary perturbation theory. For example, consider the perturbative series for the lepton 
magnetic anomalous moment: 
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a t = 



a(Q*) 



+ r 2 




+ r 3 




+ ••• 



(2.2) 



the values Q* = e~ 5 / 4 m£, etc., can be determined either by the explicit insertion of the run- 
ning coupling into the integrand of the Feynman amplitude and the mean value theorem, or 
equivalently, by simply requiring that the coefficients r n be independent of /. (Light-by-light 
scattering contributions are not related to coupling constant renormalization and thus enter 
explicitly in the order a 3 coefficient.) Thus the formula for the anomalous moment using 
the running coupling is form invariant, identical for each lepton £ = e, //, r, since the depen- 
dence on lepton vacuum polarization insertions is implicitly contained in the dependence 
of the running coupling constant. These examples are illustrations of the general principle 
that observables such as the anomalous moments can be related to other observables such 
as the heavy lepton potential V(Q 2 ) = —4ira(Q 2 )/Q 2 , which can be taken as the empirical 
definition of the on-shell scheme usually used to define ot(Q 2 ). 

The same procedure can easily be adapted |7j to non-Abelian theories such as QCD. One 
of the most useful observables in QCD is the heavy quark potential, since it can be com- 
puted in lattice gauge theory from a Wilson loop, and it can be extracted phenomenologically 
from the heavy quarkonium spectrum. If the interacting quarks have infinite mass, then all 
radiative correction are associated with the exchange diagrams, rather than the vertex cor- 
rections. It is convenient to write the heavy quark potential as V(Q 2 ) = —AhCfoiv{Q)/Q 2 - 
This defines the "effective charge" ay{Q 2 ) where by definition the "self-scale" Q 2 = —t is 
the momentum transfer squared. The subscript V indicates that the coupling is defined 
through the potential. 

In fact, any perturbatively-calculable physical quantity can be used to define an effective 
charge || by incorporating the entire radiative correction into its definition; for example 



where R° is the Born result and Q 2 = s = E 2 m is the annihilation energy squared. An 
important result is that all effective charges C(a(Q) satisfy the Gell-Mann-Low renormaliza- 
tion group equation with the same (3q and j3\\ different schemes or effective charges only 



Re+e- (Q 2 ) = Re 



(Q 2 ) 1 + 



ai?(Q) 



(2.3) 



'e+e 



TT 



S 



differ through the third and higher coefficients of the (3 function. Thus, any effective charge 
can be used as a reference running coupling constant in QCD to define the renormalization 
procedure. More generally, each effective charge or renormalization scheme, including MS, 



is a special case of the universal coupling function a(Q,/3 n ) flclfl. Peterman and Stiickelberg 



have shown JO]] that all effective charges are related to each other through a set of evolution 
equations in the scheme parameters (3 n . Physical results relating observables must of course 
be independent of the choice of any intermediate renormalization scheme. 

Let us now consider expanding any observable or effective charge ou(Qa) i n terms of 
ay ■ 



1 + (C + D VP f)^ + 

71 



(2.4) 



Since ay sums all vacuum polarization contributions by definition, no coefficient in the series 
expansion in ay can depend on /; i.e. all vacuum polarization contributions are already 
incorporated into the definition of ay. Thus we must shift the scale \x in the argument of 
ay to the scale 0] Qy = e 3Dvp ^/i: 



ola{Qa) = &v{Q\ 



* . A/V OLy 

1 + r x 1 h 

7T 



(2.5) 



where rf = C + (33/2) Dyp is the next-to- leading coefficient in the expansion of the 
observable A in scheme V. Thus the relative scale between the two observables A and V, 
^A/v = Qa/Qvi is fixed by the requirement that the coefficients in the expansion in ay 
scheme are independent of vacuum polarization corrections [f| . Alternatively, one can derive 
the same result by explicitly integrating the one loop integrals in the calculation of the 
observable A using ay(£ 2 ) in the integrand, where £ 2 is the four-momentum transferred 
squared carried by the gluon. (In practice one only needs to compute the mean-value of 
£n£ 2 = £nQ v One can eliminate the / vacuum polarization dependence that appears 

in the higher order coefficients by allowing a new scale to appear in each order of perturbation 
theory. In practice, often only the leading order commensurate scale is required in order to 
test PQCD to good precision. H 



We can compute other observables B and even effective charges such as c^Mg as an 
expansion in ay scheme: 



chb{Qb) = av(Q\ 



1 . B V OLy 

1 + r 1 h 

7T 



(2.6) 



where Qy = Qb/^b/v, and again r± must be independent of vacuum polarization contri- 
butions. We can now substitute and eliminate ay(Qy) : 



Ob{Qb) = Oi A (QA) 



* . B/A OiA , 

1 + r t ' h 

7T 



(2.7) 



where Qb/Qa = ^b/a = ^b/v/^a/v, and rf^ 4 = r^ V — r^ V . Note also the symmetry 
property Xb/a^a/b — 1- Alternatively, we can compute the commensurate scale Qa = 
Qb/^b/a directly by requiring r l to be /-independent. The result is in agreement with 
the transitivity rule: the BLM procedure for fixing the commensurate scale ratio between 
two observables is independent of the intermediate renormalization scheme. 

The scale-fixed relation between the heavy quark potential effective charge and a^jg is 
a v(Q) — a Ms ( e ~ 5 ^ 6< 3)[l — 2(aj^g/7r) + - ■ ■] 0. (The one-loop calculation oiay in MS scheme 
is given in Ref. JT5] .) 

The transitivity and symmetry properties of the commensurate scales are the scale trans- 
formations of the renormalization "group" as originally defined by Peterman and Stiickelberg 



TT| . The predicted relation between observables must be independent of the order one makes 
substitutions; i.e. the algebraic path one takes to relate the observables. It is important 
to note that the PMS method, which fixes the renormalization scale by finding the point 



of minimal sensitivity to fi, does not satisfy these group properties ||10|| . The results are 



chaotic in the sense that the final scale depends on the path of applying the PMS procedure. 
Furthermore, any method which fixes the scale in QCD must also be applicable to Abelian 
theories such as QED, since in the limit of small number of colors Nc — > the perturbative 
coefficients in QCD coincide with the perturbative coefficients of an Abelian analog of QCD 

The commensurate scale relations given in Eq. ( |2.7| ) provide a new way to test QCD: 
One can compare two observables by checking that their effective charges agree both in 
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normalization and in their scale dependence. The ratio of commensurate scales Xa/b is fixed 
uniquely: it ensures that both observables A and B pass through heavy quark thresholds at 
precisely the same physical point. Calculations are often performed most advantageously in 
MS scheme, but all reference to such theoretically-constructed schemes may be eliminated 
when comparisons are made between observables. This also avoids the problem that one need 
not expand observables in terms of couplings which have singular or ill-defined functional 
dependence. 

The physical value of the commensurate scale in ay scheme reflects the mean virtuality 
of the exchanged gluon. However, in other schemes, including MS, the argument of the 
effective charge is displaced from its physical value. The relative scale for a number of 
observables is indicated in Table I. For example, the physical scale for the branching ratio 
T — > 7X when expanded in terms of ay is (l/2.77)Mx ~ (l/3)Mx, which reflects the fact 
that the final state phase space is divided among three vector systems. (When one expands 
in MS scheme, the corresponding scale is 0.157M~f .) Similarly, the physical scale appropriate 
to the hadronic decays of the r] b is (1/1.67)M % ~ (1/2)M % . 

Table I 

Leading Order Commensurate Scale Relations 
a^(0.435Q) 



aJ1.67Q) 



a T (2.77Q) 



a T (1.36Q) ^ a v {Q) « * a*(0.614Q) 



a GLS (1.18Q) 

^a A f 2 (0.904Q) 



a sl (1.18Q) 



After scale-fixing, the ratio of hadronic to leptonic decay rates for the T has the form 
T(T -> hadrons) _ 10(tt 2 - 9) a|g(0.157M T 



r(T -> fi+fi-) 8l7re 2 b a 2 QEI) 



1 _ 14.0(5) ^MS + . 

7T 



(2.8) 
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10(vr 2 -9) c4(0.363M T ) 



8lTre 2 b a 2 QED 



1-8.0(5) — + 



(2.9) 



Thus, as is the case of positronium decay, the next to leading coefficient is very large, and 
perturbation theory is not likely to be reliable for this observable. On the other hand, 
the commensurate scales for the second moment of the non-singlet structure function M2 
and the effective charges in the Bjorken Sum Rule (and the Gross-Llewellyn-Smith Sum 
Rule) are not far from the physical value Q when expressed in ay scheme. At large n the 
commensurate scale for M n is proportional to 1/y/n, reflecting the fact that the available 
phase-space for parton emission decreases as n increases. In multiple- scale problems, the 
commensurate scale can depend on all of the physical invariants. For example, the scale 
controlling the evolution equation for the non-singlet structure function depends on xbj 
as well as Q |16 |. In the case of inclusive reactions which factorize at leading twist, each 
structure function, fragmentation function, and subprocess cross section can have its own 
scale. 

The commensurate scale relations between observables can be tested at quite low mo- 
mentum transfers, even where PQCD relationships would be expected to break down. It is 
likely that some of the higher twist contributions common to the two observables are also 
correctly represented by the commensurate scale relations. In contrast, expansions of any 
observable in cty^ (Q) must break down at low momentum transfer since a^jg (Q) becomes 
singular at Q — A^. (For example, in the 't Hooft scheme where the higher order /3 n — 
for n = 2, 3, ... , aj^(Q) has a simple pole at Q — Aj^g.) The commensurate scale relations 
allow tests of QCD without explicit reference to schemes such as MS. It is thus reasonable 
to expect that the series expansions are more convergent when one relates finite observables 
to each other. 



III. NEXT-TO-LEADING-ORDER BLM FORMULAS 



We will now generalize the BLM procedure to relate physical observables beyond leading 
order ^]. As we will see, each order in the expansion will acquire its proper commensurate 
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scale. Consider the expansion series of a physical effective charge ai(Q)/n in terms of 
another physical effective charge a 2 (Q)/n 



«i(Q) a 2 (Q) , u f\( MQ)' 



Ti 



7T 



7i 



+ (Ci2 + D 12 f + E 12 f 



'a 2 (QY 



+ 



7T 



(3.1) 



According to the BLM ansatz, we can reorganize this series and resum higher order terms 
that are induced by running of the coupling constant effects. We can perform this procedure 
order by order, and absorb these higher order terms into the renormalization scale of each 
order. Postponing the justification, we eventually obtain a series of the form 



Qi(Q) _ a 2 (Q*) ] A (^K 



C 



7i 



7i 



71 



12 



'a 2 (Q***) 



+ 



7T 



(3.2) 



Where the running coupling constant effects have been resummed into the renormalization 
scales Q*, Q** and Q***. 

We now analyze the steps involved in this process of resummation. First of all, let us 
make an important observation: the Bi 2 f term and the E\ 2 f 2 term come exclusively from 
the one-particle-irreducible vacuum polarization graphs with one and two fermion loops. 
These fermionic contributions belong to the running coupling constant effects, and should 
be fully absorbed into the renormalization scale Q*. 

There are also some gluonic contributions in the running coupling constant effects, and 
consequently some part of A 12 ,C\ 2 and D X2 should also be absorbed into Q*. The exact 
amount is dictated by the behavior of the running coupling constant. The running of a 2 (/j.) 
for a general SU(N) group can be characterized by 



a 2 (jj) a 2 (n 



o 



1 /H 



71 



7Y 
1 

16 



C 



4 V 3 
11 I 



A 



- ! i -C A - -77 ) £n 2 (^ 



'a 2 (no) 



71 



A, 



34 



C 



20 



C A + C F ) Tf 



in 



' ot 2 {n G 



71 



+ 



(3.3) 
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In the above formula, Ca = N, Cf = {N 2 — 1)/2N are the quadratic Casimirs of the adjoint 
and the fundamental representations. T is the normalization of the trace of generators of 
the fundamental representation: Tr(t a t b ) = T5 ab . Conventionally T is chosen to be 1/2. 

Let us use a 2 (fi) as the strong coupling constant in the calculation of the physical effective 
charge «i(Q). If we identify the renormalization scale /i with the momentum transfer of 
the exchanged gluons in the leading order (LO) Feynman diagrams, and set /xo = Q, the 
contribution from these diagrams will have the form 



«i(Q) 



7T 



a 2 (Q) i a i CA _^ Tf \ Li f^(Q)' 



lo 71 4V3~~ 3"V" X V 71 ) 
1 f /ll 4 \ 2 

34 . /20_ , \ \ 1 fa 2 (Q)' 1 



T CX- [ Y C A + C F) Tf )Ll ^ + 



(3.4) 

The quantities L x and L 2 are two numerical constants, which may be interpreted as the 
mean- values of the £n(/i 2 /Q 2 ) and £n 2 (/i 2 /Q 2 ) terms. In general L 2 ^ L\ due to loop 
integration or due to the presence of more than one Feynman diagram involving a tree- 
level gluon propagator. This is the structure of terms that we should absorb in Eq. ( |3.1| ). 
In other words, we can shift the renormalization scale in Eq. (|3.1| ) until we fully absorb 
those higher order terms with the structure described in Eq. ( p.4|) . To make this step 
more clear, let us focus on the NLO coefficient of Eq. ( |3.1|) , which is given by A\ 2 + B\ 2 f '. 
When we perform a scale shift Q — ► Q*, we eliminate the Bi 2 f term completely, but at the 
same time, we also modify the A\ 2 term, because the net change to the NLO coefficient is 
proportional to yC^ — as indicated by the structure of the NLO coefficient in Eq. 

( p.4| ). The value of L\ is determined by the condition of eliminating the B i2 f term. Now, 
proceeding to the NNLO coefficient of Eq. ( |3.1| ), which is given by C\ 2 + D i2 f + E 12 f 2 , we 
can determined the correct value of L 2 to eliminate the E\ 2 f 2 term. The modification to 
the C\ 2 and the D\ 2 f term is given by the form of the NNLO coefficient in Eq. ( ^.4| ). The 
exact form of the scale shift Q — > Q* is given later in Eq. (|3.15 ). We should clarify that 
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the scale setting happens automatically when we use the running coupling constant inside 
the Feynman diagram calculation (like the skeleton expansion in QED). The mathematical 
process here of eliminating the B\ 2 f and terms allows us to determine the exact amount 
of the scale shift. 

After this first step, our series will look like 

?m. = ?m + Al2 ( a -m)\ {c - 12 + D[2f) (^m) 3 + . . . . ( 3. 5) 

The next step to follow is now clear. We should absorb the D' 12 f term into a new renor- 
malization scale Q** for the NLO term, since the / dependence comes from the vacuum 
polarization corrections to the NLO term. Naturally part of the C[ 2 term should also be 
absorbed. From the form of the running coupling constant in Eq. (|3.3|) , we know that the 
absorbed term should be proportional to yCU — After this procedure we arrive to the 
form indicated in Eq. ( |3.2| ). To the order considered here, we do not have enough informa- 
tion to set the scale Q*** of the next-to-next-to-leading order (NNLO) coupling constant. A 
sensible choice is Q*** = Q**, since this is the renormalization scale after shifting the scales 
in the second step of the BLM procedure. 

In practice, most physical observable in perturbative QCD are computed in the MS 
scheme, with the running coupling constant fixed at the physical scale of the process. Specif- 
ically, if the perturbative series for ai(Q)/n and a2(Q)/n are 

<*i(Q) = a Ms(Q) + , Ai + Bi fs ( a us(Q) Y 

7T 7T \ TT J 

+{d + DJ + E,f) (^P) 3 + • • • , (3.6) 



7T 7T \ TT J 

+ (C 2 + D 2 f + E 2 f) (^P) 3 + " " " , (3.7) 
then the coefficients A 12 , B 12 , Ci 2 , D 12 and E x % of Eq. ( |3.1| ) are given by 
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A 12 = Ai- A 2 , 
B12 — B\ — B 2 , 

Cia = Ci - C 2 - 2(Ai - A 2 )A 2 , 

D 12 = D X -D 2 - 2(A 1 B 2 + A 2j Bx) + 4A 2 5 2 , 

E 12 = Ei — E 2 — 2{Bi — B 2 )B 2 , 

and the NLO BLM formulas (for a general SU(N) group) are given by 



(3.8) 
(3.9) 
(3.10) 
(3.11) 
(3.12) 



A 12 
C\ 2 
Q* = Q exp 



A\ 2 + 



Q** = Qexpl ^^12 



Bio 

2T 



-A' 1 



9 



8T 2 V 3 
1 



^C A - t T f) (B\ 2 - E 



12 



a 2 (Q) 



7T 



-{hC A + 3C F )B 12 + D 12 + YY~ E ^ 



(3.13) 
(3.14) 
(3.15) 

(3.16) 



Notice the presence of a 2 {Q)/rt in the expression of Q*. In general Q* will itself be a 
perturbative series in a 2 (Q)/n. This fact has first been pointed out in Ref. 0], and also 
been further explored by Grunberg and Kataev on the extension of the BLM approach [|17| . 
We have exponentiated the perturbative series, since physically the renormalization scale Q* 
should always be positive. To the order considered here, the scale for the coupling constant 
in Eq. (|3.15| ) is not well-defined, but can be chosen to be Q. This intrinsic uncertainty 
is similar to the Q*** scale uncertainty of the NNLO term in Eq. (|3.2|) , and can only be 
resolved by going to the next-higher order. 

For a general gauge field theory, it is interesting to point out that the scale setting 
procedure described here leads to the correct expansion series coefficients in the "conformal 
limit" HH]. The conformal limit is defined by /3o, Pi — » 0, and can be reached, for instance, by 
adding enough spin-half and scalar quarks. Since all the running coupling effects have been 
absorbed into the renormalization scales, the scale setting method described here correctly 
reproduces the expansion coefficients in this limit. It should be pointed out that other scale 
setting methods in general do not guarantee this feature. 
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The application of the scale setting formulas to currently available NNLO QCD quantities 
gives some very interesting results. The following is a list of some effective charges known 
to the NNLO. 

1) olr(Q)/it: the effective charge obtained in total hadronic cross section in e + e~ annihila- 



tion, defined by [flj 



R(Q)=3Y,Q 2 f 

f 



1 + 



<xr(Q) 



71 



(3.17) 



The perturbative series of o>r(Q)/7t is (using T = 1/2 for the trace normalization) 
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(3.18) 



The term containing (J2 / Q f) 2 1 / Q} arises from light-by-light diagrams. The dimension 
of the quark representation is d(R), which usually is N for SU(N). For QCD we have 
d abc d abc = 40/3. It is interesting to point out that the vector part of the Z hadronic decay 
width [T3|] shares the same effective charge as R{Q). 

2) a T (M T )/ir: the effective charge obtained from the perturbative hadronic decay rate of the 



r lepton, defined by |l| 



R T = 3 



1 + 



a T (M T ) 



71 



(3.19) 



The perturbative series of a T (M T )/n is 

<*AQ) = <WQ) + ( a m(Q) \ 2 
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(3.20) 



The appropriate number of flavors to be used here is / = 3. There is no light-by-light 
contributions in a T /ir. 

3) oi gi (Q)/iv: the effective charge obtained from the Bjorken sum rule for polarized electro- 



production, defined by ]K 



dx 



g?{x,Q 2 )-g{ n {x,Q 2 
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(3.21) 



Notice that there are different normalization conventions in the literature for the polarized 
structure functions. Here we follow the normalization given in Ref . |19| . For a recent review 



on this sum rule, see Ref. | 20| . The perturbative series of a gi (Q)/ir is 

2 r: 



Qgi (Q) = a Ms(Q) + ( a Ms(Q) 



7i 



71 



71 



2 *c A - 7 -c F - l -f 

12 8 3\ 



+ 



a- 



MS 



7T 



V 648 



1241 11 



432 9 



+ ^fCs ) c A c F + -^cl 



1 

32 



+ 



115 

648 



3535 1 5 \ /133 5 

"1296 " 2 Cs + 9 Cs J ° A + l864 + 18 Cs ' Cl 



f 



(3.22) 



For a 9l /7T we do not have light-by-light contributions, either. 

4) ap 3 (Q) /tt: the effective charge obtained from the Gross-Llewellyn Smith sum rule, defined 



by 



dx 



Fz p (x, Q 2 ) + Fs p (x, Q 2 



= 6 


x <*F Z {Q) 




71 



(3.23) 
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As pointed out in Ref. |IJ, the effective charge OF 3 /7r differ from a gi /7i only by the light- 
by-light contributions. The perturbative series of OLp 3 j'K is given by the perturbative series 
of a 9l /vr plus the additional term in the following formula. 



a F 3 (Q) _ a gi {Q) ( a m {Q) 



3 d abc d abc 



77 



77 



77 



11 + ks I /■ 



(3.24) 



/ C F N C V 144 ' 6' 
For QCD we have Cp = 4/3, the number of colors is Nq = 3, and d abc d abc = 40/3. Since the 
radiative corrections to the Bjorken sum rule are identical to those of the Gross-Llewellyn- 
Smith sum rule — up to small corrections of order a^(Q 2 ), a basic test of QCD can be made 
by considering the ratio of the Gross-Llewellyn-Smith and Bjorken sum rules: 

1 ■£ dx \f?{x, Q 2 ) + F? (x, Q 2 ) 
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Since the Regge behavior of the two sum rules is similar, the empirical extrapolation to 
e — > should be relatively free of systematic error. Moreover, PQCD predicts 



RgllsibAQ\ e 0) = 1 + O (a 3 s (Q)) + O 



'A 



QCD 



Q' 2 



(3.26) 



i.e., hard relativistic corrections to the ratio of the sum rules only enter at three loops. Thus 
measurements of the ratio of the sum rules could provide a remarkably complication-free 
test of QCD - any significant deviation from RcLLS/BjiQ 2 , e — > 0) = 1 must be due to higher 
twist effects which should vanish rapidly with increasing Q 2 . 

5) otF 1 (Q)/7T: the effective charge obtained from the Bjorken sum rule for deep- inelastic 



neutrino-nucleon scattering, defined by JH] 



dx 
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(3.27) 



The perturbative series of cxf 1 (Q)/'k is 
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For aF 1 /ir we do not have light-by- light contributions, either. 

As the first example of a beyond leading order commensurate scale relation, we shall 
express a T (M T )/iT in terms of a R (Q)/ir. The appropriate number of flavors in this case is 
/ = 3, because r decay occurs below the charm threshold. [Incidentally, the light-by-light 
contribution in a R {Q)/n vanishes for the three flavor case.] The application of the NLO 
BLM formulas leads to the following commensurate scale relation 



c*t(M t ) = a R {Q*) 

TT TT 

Q* = M T exp 



19 169 a R (M T ) 
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(3.29) 
(3.30) 



24 m 

Notice that all the (3,(5 an d tt 2 terms present in the perturbative series of a R (Q)/TT and 
a T (M T )/iT have disappeared when we related these two physical observables directly. Notice 
also the vanishing NLO and NNLO coefficient in Eq. ( |3.29| ). That is, up to the NNLO, the 
two effective charge are simply related by a BLM scale shift. 

As the next example, let us express a gi {Q)/iT in terms of oi r {Q)/tt. We will leave the / 
dependence explicit. The application of the NLO BLM formulas leads to 



a 9l (Q) _ a R (Q* 
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As explained previously, the scale Q*** in the above expression can be chosen to be Q** . 
Notice that aside from the light-by-light contributions, all the £3, (5 and tt 2 dependencies 
have been absorbed into the renormalization scales Q* and Q** . Understandably, the tt 2 term 
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should be absorbed into renormalization scale since it comes from the analytical continuation 
of QCD coupling constant from the spacelike region to the time-like region. However, at 
present we do not have an full understanding for the disappearance of the £3 and £5 terms. 
For the three flavor case, or neglecting the light-by-light contribution, the series simplifies 

to 

a gi (Q) = &r(Q*) ~ °?r{Q**) + &r(Q~), (3-34) 



where 



3C 

a gi (Q) = -^a gi (Q), (3.35) 
3C 

a R (Q) = ~^-oc R {Q). (3.36) 



Observe the asthonishing simplicity of the expansion series in Eq. (|3.34j) . These last 
formulas suggest that for the general SU(N) group the natural expansion parameter is a. 
The use of a also makes explicit that the same formula is valid for QCD and QED. That 
is, in the limit No — > the perturbative coefficients in QCD coincide with the perturbative 
coefficients of an Abelian analog of QCD ||. 

Broadhurst and Kataev have recently observed a number of interesting relations between 
ozr(Q) and a gi (Q) (the "Seven Wonders") J22|. Our simple result here reinforces the idea 
of a "secret symmetry" between «_r(Q) and a gi (Q). We see that this hidden simplicity is 
only revealed after the application of the NLO BLM scale setting procedure. 

The application of the NLO BLM formulas to related the other effective charges presented 
here give results with similar simplicity (see Appendix), except for those cases involving a Fl . 
For instance, the obtained relation for a T (M T ) and ap 1 is 
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(3.38) 
(3.39) 



2 15^ 3 

We see that the it 2 terms is absorbed into the renormalization scale Q*, but the resulting 
coefficients in the expansion series are not as simple as in the previous examples. In partic- 
ular, we observe that the (3 and £5 terms persist in the NNLO coefficient, which also has a 
Ca dependence. 

When the scale setting formulas are applied to relate physical effective charges to the 
somewhat artificial MS scheme coupling constant, no comparable simplicity is observed. For 
instance, the application of the formulas to relate &r{Q) to ot^g(Q) leads to 
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These last expressions clearly do not exhibit the simplicity of the cases shown in Eq. ( |3.29| ) 
or Eq. (pi). 



IV. CONCLUSION 



The problem of the scale ambiguity of PQCD predictions has plagued attempts to make 
reliable and precise tests of the theory. In this paper we have shown how this problem 
can be avoided by focussing on relations between experimentally-measurable observables. 
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The conventional MS renormalization scheme serves simply as an intermediary between 
observables. For example, consider the entire radiative corrections to the annihilation cross 
section expressed as the effective charge a R (Q) where Q = y/s: 

ocr(Q) 



R(Q) = 3J2Q} 
f 



1 + 



TT 



(4.1) 



Similarly, we can define the entire radiative correction to the Bjorken sum rule as the effective 
charge a gi (Q) where Q is the lepton momentum transfer: 



dx 



g?(x,Q 2 )-gr(x,Q 2 



9a 



9v 



1 



a gi (Q) 



TT 



(4.2) 



We now use the known expressions to three loops in MS scheme and choose the scales 
Q* and Q** to re-sum all quark and gluon vacuum polarization corrections into the running 
couplings. The value of these scales are the physical values which insure that each observable 
passes through the heavy quark thresholds at their respective commensurate physical scales. 
The final result is remarkably simple: 



a gi (Q) _ a R (Q* 



+ 



+ 



(4.3) 



TT TT \ TT J \ TT J 

A fundamental test of QCD is to verify empirically that the observables track in both 
normalization and shape as given by these relations. The coefficients in the series (aside 
for a factor of Cp which can be absorbed in the definition of a s ) are actually independent 
of color and are the same in Abelian, non-Abelian, and conformal gauge theory. The non- 
Abelian structure of the theory is reflected in the scales Q* and Q**. The commensurate scale 
relations thus provide fundamental tests of QCD which can be made increasingly precise and 
independent of any scheme or other theoretical convention. 

We have also presented in this paper a number of other commensurate scale relations 
using the extension of the BLM method to the next-to-leading order. We have shown that 
in each case the application of the NLO BLM formulas to relate known physical observables 
in QCD leads to results with surprising elegance and simplicity. 

In principle, commensurate scale relations allow tests of perturbative QCD with higher 
and higher precision as the perturbative expansion grows. They also provide a new way 
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to specify QCD phenomenology. Because they relate observables, the commensurate scale 
relations are convention-independent; i.e., independent of the normalization conventions 
used to define the color SU (N) matrices, etc. Since the ambiguities due to scale and scheme 
choice have been eliminated, one can ask fundamental questions concerning the nature of 
the QCD perturbative expansions, e.g., whether the series is convergent or asymptotic due 
to renormalons, etc. ||23|| . The precision of the leading twist series will also allow sensitive 



tests for higher twist contributions to physical observables. 

We emphasize that any consistent renormalization scheme, MS, MS, MOM, etc., with any 
arbitrary choice of renormalization scale \x, can be used in the intermediate stages of analysis. 
The final result, the commensurate scale relation between observables, is guaranteed to be 
independent of the choice of intermediate renormalization scheme since the BLM procedure 
satisfies the generalized renormalization group properties of Peterman and Stuckelberg. An 
important computational advantage is that one only needs to compute the /-dependence of 
the higher order terms in order to specify the lower order scales in the commensurate scale 
relations. In many cases, the series coefficients in the commensurate scale relations can be 
determined from the corresponding Abelian theory; i.e. Nq — > 0. 

The BLM method and the commensurate scale relations presented here can be applied 
to the whole range of QCD and standard model processes, making the tests of theory much 
more sensitive. The method should also improve precision tests of electroweak, supersym- 
metry and other non-Abelian theories. One of the most interesting and important areas of 
application of commensurate scale relations will be to the hadronic corrections to exclusive 
and inclusive weak decays of heavy quark systems, since the scale ambiguity in the QCD ra- 
diative corrections is at present often the largest component in the theoretical error entering 
electroweak phenomenology. 

The commensurate scale relations for some of the observables discussed in this paper 
(or, a T , a gi and of 3 ) are universal in the sense that the coefficients of a s are independent of 
color; in fact, they are the same as those for Abelian gauge theory. Thus much information 
on the structure of the non-Abelian commensurate scale relations can be obtained from much 
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simpler Abelian analogs. In fact, in the examples we have discussed here, the non-Abelian 
nature of gauge theory is reflected in the /3-function coefficients and the choice of second- 
order scale Q**. The lack of convergence of the non-Abelian theory such as renormalon 
behavior could show up as a progressive decrease of the higher order commensurate scales. 
The coefficients in these commensurate scale relations are simply ±1 and and suggest that 
the underlying relation between observables without light-by-light contributions is possibly 
a geometric series. Note that the relative correction due to the three loop corrections is 
±(a s /7r) 2 or less than 1% for a s < 0.3. 

A natural procedure for developing a precision QCD phenomenology is to choose one 
effective charge as the canonical definition of the QCD coupling, and then predict all other 
observables in terms of this canonical measure. Ideally, the heavy quark effective charge 
a v{Q 2 ) could serve this central role since it can be determined from both the quarkonium 
spectrum and from lattice gauge theory. However, in order for this effective charge to be 
useful in practice, it will be necessary to compute the relation of the heavy quark potential 
to other schemes through three loops. At present, the most precisely theoretically and 
empirically known effective coupling is a#(Q 2 ), as determined from the annihilation cross 
section; thus it is natural to use it as the standard definition. 

Alternatively, one can follow historical convention and continue to use the MS scheme 
as an intermediary between observables. For definiteness, we can define the MS scheme 
as having all /3 n = beyond n = 1. The commensurate scale relations such as Eq. ( p.l5|) . 



Q3.16Q then unambiguously specify all of the scales Q*, Q**, etc. required to relate to the 



observables. The intrinsic QCD scale will then be unambiguously encoded as Aj^g. However, 
there is am intrinsic disadvantage in using a^(Q) as an expansion parameter: the function 
a Ms(Q) nas a simple pole at Q = A^jg, whereas observables are by definition finite. 



The BLM scale has also recently been used by Lepage and Mackenzie |nj and their 
co-workers to improve lattice perturbation theory. By using the BLM method one can elim- 
inate lattice in favor of ay thus avoiding an expansion with artificially large coefficients. 
The lattice determination, together with the empirical constraints from the heavy quarko- 
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nium spectra, promises to provide a well-determined effective charge ay{Q) which could be 
adopted as the QCD standard coupling. 

After one fixes the renormalization scale fi to the BLM value, it is still useful to com- 
pute the logarithmic derivative of the truncated perturbative prediction din p^jdin\x at the 
BLM-determined scale. If this derivative is large, or equivalently, if the BLM and PMS 
scales strongly differ, then one knows that the truncated perturbative expansion cannot be 
numerically reliable, since the entire series is independent of //. Note that this is a necessary 
condition for a reliable series, not a sufficient one, as evidenced by the large coefficients in the 
positronium and quarkonium decay widths which appear when the scales are set correctly. 
In the case of the two and three jet decay fractions in e + e~ annihilation, the BLM and PMS 
scales strongly differ at low values of the jet discriminant y. Thus, by using this criterion, 
we establish that the perturbation theory must fail in the small y regime, requiring careful 
resummation of the a s £ny series. (A more detailed discussion of the sensitivity of the jet 
fractions to scale choice and jet clustering schemes is given in Ref. 0.) 

However, if we restrict the analysis to jets with invariant mass M. < y/ys, with 0.14 > 
y > 0.05, then we have an ideal situation, since both the PMS and FAC scales nearly 
coincide with the BLM scale when one computes jet ratios in the MS scheme (See Fig. 1.) 
i.e., the renormalization scale dependence in this case is minimal at the BLM scale, and the 
computed NLO (next-to-leading order) coefficient is nearly zero. In fact, Kramer and Lampe 
H find that the BLM scale and the NLO PQCD predictions give a consistent description of 
the LEP 2-jet and 3-jet data for 0.14 > y > 0.05 at the Z. Neglecting possible uncertainties 
due to hadronization effects, this allows a determination of a s with remarkably small error: 
g a m (M z ) = 0.107 ± 0.003, which corresponds to = 100 ± 20 MeV. 

The central principle we have used in our analysis is the fact that vacuum polarization 
contributions are summed by the running coupling constant in gauge theory. The argument 
of the running coupling constant is then fixed by the requirement that all fermion vacuum 
polarization is resummed into a s , rather than appear in the coefficients. The fact that NLO 
correction to the scale Q* is proportional to fii is consistent with the Peterman-Stuckelberg 
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renormalization group analysis and is crucial for applying this method to higher order. 
We have also seen that this scale-setting procedure leads to correct expansion coefficients 
in the conformal limit, since the beta-function dependence has been resummed into the 
renormalization scales. The same procedure can be applied to multi-scale problems; in 
general, the commensurate scales Q*,Q**, etc. will depend on all of the available scales. 
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In this Appendix we present the results of applying the next-to-leading order BLM pro- 
cedure to relate the five known effective charges aR(Q),a T (Q),a gi (Q) aF 3 (Q) and ofi(<5)- 
Notice that in principle a T (Q) is meaningful only for the three-flavor case and for Q = M T . 
However, here we will use its analytical expression and leave Q and / dependence explicit. 
We will leave the scale Q*** in the following formulas unspecified. This scale in general can 
be chosen to be Q**, or in the absence of the NLO term, Q*. We use Tr (t a t b ) = \5 ab for the 
normalization of the trace generators in this Appendix. 
or(Q) in terms of a T (Q). 
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renormalization scale as given by A^/Ay leads to a vanishing NLO coefficient, which 
is equivalent to the FAC criterion. In contrast, in the BLM approach, only the / part 
vanishes. 
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Celmaster and Stevenson |2(| have observed that the NLO coefficient r^ A in an expansion 



series such as Eq. Q2.5|) depends on the quantity A used in the right-hand-side, and 
they have expressed concern that this fact implies a "scheme dependence" of the BLM 
procedure. This is obviously not a defect: when one expands an observable B in terms of 
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